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Abstract 

I perform an independent analysis of radio Doppler tracking data from the Pioneer 10 spacecraft 
for the time period 1987-1994. All of the tracking data were taken from public archive sources, and 
the analysis tools were developed independently by myself. I confirm that an apparent anomalous 
acceleration is acting on the Pioneer 10 spacecraft, which is not accounted for by present physical 
models of spacecraft navigation. My best fit value for the acceleration, including corrections for 
systematic biases and uncertainties, is (8.60 dt 1.34) x 10~ 8 cm s -2 , directed towards the Sun. 
This value compares favorably to previous results. I examine the robustness of my result to 
various perturbations of the analysis method, and find agreement to within ±5%. The anomalous 
acceleration is reasonably constant with time, with a characteristic variation time scale of > 70 yr. 
Such a variation timescale is still too short to rule out on-board thermal radiation effects, based 
on this particular Pioneer 10 data set. 
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I. INTRODUCTION 


Measurements of spacecraft motions in the solar system can be used as tests of gravitation 
and relativity. Recently, Anderson et al. [l, 2] have presented the discovery of an anomalous 
effect seen in radio tracking data from the Pioneer 10 spacecraft. When interpreted as a 
Doppler shift, this anomalous effect corresponds to a constant acceleration, directed towards 
the Sun, of approximately (8 ± 1) x 10 -8 cm s -1 . Anderson et al. (hereafter A02) found 
that the anomalous effect could not be explained by previously known physics or spacecraft 
properties. 

The discovery of the anomaly has stimulated numerous efforts to explain it. Some of the 
explanations involve “new physics,” such as modified gravity or dark matter, while other 
explanations invoke a change in the physical properties of the Pioneer spacecraft, such as a 
asymmetric radiation profile. I have considered a third avenue of exploration, which is to 
test the analysis procedure for flaws. In this paper, I present an independent analysis of the 
Pioneer 10 trajectory and search for an anomalous acceleration. 

A02 studied radio tracking data from four deep space missions: Pioneer 10, Pioneer 11, 
Ulysses, and Galileo. All four of these missions showed suggestions of an anomalous accel- 
eration of order 10~ 7 cm s~ 2 . However, A02 considered the determination of the anomalous 
acceleration of the Pioneer 10 spacecraft to be the most reliable. Therefore I have focussed 
exclusively on the Pioneer 10 data for my analysis. Once I had verified the presence of an 
anomalous acceleration, I tested the result for robustness in several different ways. 

All of the procedures discussed in this paper were developed by myself [9] and written 
in the Interactive Data Language (IDL) [7, 8]. During the development I had only minimal 
contact with the A02 group authors, and as I detail below, these contacts had a minimal 
impact. Thus, I consider this work to be an independent test of the analysis by A02. I have 
analyzed a subset of the Pioneer 10 tracking data that is available from the public archives, 
which is, of course, the same data that A02 used in their analysis. The time coverage of my 
analysis (years 1987-1994) is most comparable to that of the original discovery presented 
in A02. By necessity, many of the analysis procedures I developed will be at least similar 
to those of A02, but I attempt to extend the analysis by considering additional models, 
including spacecraft spin, maneuvers, and time- variation of the anomalous acceleration. 

The contents of the paper are as follows. Section II briefly describes the Pioneer 10 
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- spacecraft and the Deep Space Network systems. Section III presents the methods I used 
to acquire the data and perform initial filtering, while Sec. IV describes the analysis and 
modeling techniques that I employed. The results of the tracking and uncertainty analyses 
are presented in Sections V and VI. This is followed by a short discussion and conclusion 
in Sections VII and VIII. 

II. SPACECRAFT AND COMMUNICATIONS 

I provide only a cursory description of the Pioneer spacecraft, and Pioneer and DSN 
communications systems here. I refer the reader to A02, and references therein, for a more 
complete description. 

A. Pioneer 10 Spacecraft 

Pioneer 10 was launched on 2 March 1972 and, after an encounter with Jupiter, has 
followed a hyperbolic escape trajectory from the solar system. On 1 January 1987, the 
spacecraft was approximately 40 A.U. from the Sun, and receding with a nearly constant 
velocity of 12.8 km s -1 . 

The main physical features of the Pioneer 10 spacecraft are the parabolic high gain 
antenna, with a radius of 137 cm; the instrument compartment, which faces the direction of 
travel; and the two radioisotope thermoelectric generators (RTGs), which are attached to 
the compartment by booms of 300 cm length. The Pioneer 10 spacecraft is spin stabilized, 
and spins at a rate of ~ 5 rpm. The spin axis is aligned with the high gain antenna axis, 
which is designed to point towards the Earth, opposite the direction of travel. 

The transponder on board Pioneer 10 functions in the S-band. The uplink signal from 
Earth is received by the high gain antenna at 2.11 GHz, while the downlink signal to Earth 
is transmitted at a frequency close to u 0 = 2.292 GHz. I have used only data where the 
spacecraft communications system was operated in “coherent” mode. In this mode, the 
spacecraft retransmits a downlink signal, phase-locked to the uplink, with an exact frequency 
turnaround ratio of 240/221. 
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B. DSN Communications 

Communications with Pioneer 10 are accomplished using the Deep Space Network (DSN). 
The DSN maintains antenna complexes at Canberra, Australia; Goldstone, California; and 
Madrid, Spain. Radio tracking observations are normally obtained in the course of mission 
operations for the purposes of spacecraft navigation. The two basic types of tracking are 
ranging, which directly measures the spacecraft distance via round trip light travel time; 
and Doppler, which measures indirectly the range rate, or relative velocity along the line of 
sight. Tracking passes were obtained more or less regularly spaced throughout the time range 
1987-1994, and within individual years, although there were apparently some campaigns of 
more intensive tracking coverage. 

During uplinks, a digitally controlled oscillator (DCO) is programmed to a precise fre- 
quency, which then drives the Exciter Assembly, whose signal is sent to the transmitting 
antenna. The uplink frequency is typically Doppler-compensated so that the frequency re- 
ceived by the spacecraft is near 2.11 GHz [3]. On the downlink leg, after being received at 
Earth, the Pioneer 10 Doppler signal is down-converted to a 1 MHz intermediate frequency. 
The Metric Data Assembly is used to accumulate a continuous count of Doppler cycles 
at this intermediate frequency, at fixed integration intervals. The Doppler count so com- 
puted represents the integrated range rate (i.e. the line of sight change in distance between 
the antenna and spacecraft). These quantities are typically differenced and, after further 
straightforward manipulations, produce the mean Doppler frequency over the integration 
interval [2, 4, 5]. 

Since the Pioneer spacecraft was 40-69 A.U. from the Earth for the period of this analysis, 
the round trip light travel time ranges from 11 to 16 hours. Thus, while it is possible for the 
same station to be used for both the uplink and downlink legs of the transmission to and 
from the spacecraft (known as “two-way” Doppler), it is more common for one station to 
transmit the uplink leg and a separate station to receive the downlink leg some hours later, 
which is known as “three-way” Doppler. As far as the analysis is concerned, both of these 
kinds of Doppler tracking data are identical and can be handled the same way. 

For spacecraft ranging, a unique repeating ranging code is modulated onto the 2 GHz 
carrier wave. Upon return from the spacecraft, the received ranging code is correlated with 
the transmitted one, and a range time delay can be computed, modulo the period of the 
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ranging code pattern. No reliable range data were available for Pioneer 10, and so I analyzed 
only the Doppler tracking data. 


III. DATA PREPARATION 

I obtained Pioneer 10 Doppler tracking data from the publicly accessible NSSDC archive. 
The JPL Radio Science group has submitted a substantial portion of the tracking data to 
the archive in the form of digital tapes which must be staged manually onto a computer 
disk. I requested data covering the time period 1987-1994. Data from beyond April 1994 is 
not available from NSSDC. In addition, there are a number of gaps in the data. The largest 
gap is due to an unreadable archive tape which covered the interval June 1990 to June 1991. 
The data are stored in a standard Archival Tracking Data File (ATDF) format [6]. 

An initial level of filtering and processing was applied to the raw ATDF records. A large 
number of records were in “one-way” Doppler mode (i.e., transmissions originating from the 
spacecraft) and were simply discarded. The integration time of the records was variable, and 
ranged from 0.1 s (so-called “high rate” Doppler), to ~ 100 s or more. The high rate data 
in particular contained a large number of samples, and to prevent over-weighting of those 
segments, I chose to accumulate the Doppler counts to intervals of at least 60 s in duration. 
I also eliminated discontinuous or noisy data, which occur preferentially at the beginnings 
and ends of tracking passes, or during noisy passes. A basic sliding 10-sample median filter 
was applied, and points more than 100 Hz from the median were discarded. 


IV. ANALYSIS OVERVIEW 

The data consist of a time series of observed frequencies at designated DSN antennae. The 
dominant variations observed in the data are the annual signature of the Earth’s motion 
in the solar system and the diurnal signature of the Earth’s rotation. The diurnal term 
contains both the motions of the receiving and the transmitting antennae involved in the 
tracking pass. Finally, of course, there is the Doppler variation introduced by the Pioneer 
10 spacecraft itself, which is the sought-after signal. 

In this section, I will follow the terminology established by A02 in identifying epoch of 
transmission and reception. The epoch of transmission from the Earth is ti, the epoch of 


5 



interaction of the signal with the Pioneer 10 spacecraft is t 2 , and the epoch of reception back 
at the Earth is t 3 . All of these times are referred to the Barycentric Dynamical Timescale 
(TDB), which is a coordinate time at the solar system barycenter. TDB is also the effective 
argument of the JPL planetary ephemerides. The 3- vectors ri, r 2 , and r 3 represent the 
positions of the corresponding antenna at the corresponding epoch, and Vi, v 2 , and v 3 
represent the velocities. The vector difference, r 12 , is defined as rq — r 2 . These vector 
quantities are measured in the solar system barycenter frame. 

The original station times in the ATDF records are referred to Coordinated Universal 
Time (UTC). When computing Earth rotation and orientation quantities, the Terrestrial 
Dynamical Time (TDT) timescale is used. Conversion between the UTC, TDT and TDB 
timescales is straightforward using standard practices [10, 11]. 

The expected frequency at the receiver at time f 3 can be expressed as 

fie — [(240/221) — 77 /spin] 7^23 + Az2 path , (1) 

where f\ is the uplink frequency at time t x as measured at the transmitter, / spin is the 
spacecraft spin frequency at time t 2 , and the ratio 240/221 is the spacecraft transponder 
turnaround ratio (note that 77 = 1 + 240/221). The factors d 12 and cf 23 embody the Doppler 
shifts of the moving spacecraft and earthbound antennae. The frequency multiplier on the 
uplink leg is 

w _ \/l-lvil 2 /c 2 (1 - f 12 • v 2 /c 2 ) ^ 

fll2 (1 - f 12 • Vl /c 2 ) Vl - H 2 /c 2 ’ 1 ' 

where the first fraction represents the relativistic Doppler shift due to the Earth motion, and 
the second due to the spacecraft motion. The unit vector r 12 points from the transmitting 
station to the spacecraft, i.e., r i2 = r 12 /r 12 . The downlink factor, d 23 , is constructed in the 
same fashion, by substituting 1—^2 and 2 — * 3. 

The final term in equation 1, Az7 path , represents additional Doppler effects caused by 
small effective path length changes, aside from those due to geometric antenna motions. 
Generally speaking, this term can be written as A^ pat h = —2 dl/dt v 0 jc, where dl/dt is the 
time rate of change of effective photon trajectory path length along the line of sight. The 
factor of 2 comes from the two legs of the round trip path. 

In this paper I consider the effective path length due to the “Shapiro” delay [12]. The 
Shapiro delay reflects the extra proper distance traveled by a photon, beyond the classical 
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geometric distance, in the Sun’s gravitational potential, as predicted by general relativity, 


4hap — (1 + 7) 


G Mr: 


In 


r x + r 2 + r i2 
n+r 2 - r 12 J ’ 


( 3 ) 


where 7 is a parameter of the parameterized post-Newtonian formulation of gravity [13, 14]. 
For general relativity, 7 = 1. On an annual timescale, the impact parameter of the photon 
trajectory increases and decreases, with a minimum distance of about 8x 10 6 km. Conversion 
to a Doppler shift is achieved by numerically differentiating equation (3), which yields an 
annual signal with amplitude ±150 mHz. As discussed further below, I do not model the 
effects of the solar corona. 

The known quantities are the receiver quantities £3, fz, and the station identification for 
each Doppler sample. In order to compute the expected frequency at the same epoch, all 
of r{x )2) 3} and V{i )2> 3j must be determined. This is especially important because even the 
time of transmission, t\, and hence the transmission frequency, /1, are not known a priori. 
Starting from the reception epoch, the spacecraft epoch t 2 is determined by solving the light 
travel time equation [8], tz — t 2 = r 2 z/c + (£ S hap)23 via an iterative process, using the known 
trajectory and rotation properties of the Earth, and a trial trajectory of the spacecraft. In 
the same way, the transmission epoch ti can be determined. The ATDF data contains a 
special record of the transmitter configurations, including the frequency of the antenna’s 
DCO, from which the transmitted frequency fi can be determined. 

The motions of the Earth center are interpolated from the JPL DE405 planetary 
ephemeris [8, 22], which is referred to the axes of the International Celestial Reference Frame 
(ICRF). The position and velocity of Earth stations with respect to the geocenter, referred 
to an inertial coordinate system such as the ICRF, must take into account Earth rotation 
and the changes in Earth orientation parameters. The apparent sidereal time, which is the 
hour angle of the Earth referred to an inertial system, is taken from Aoki et al. [10, 16]. 
Earth precession and nutation describe the motion of the Earth rotation axis with respect 
to the celestial sphere. These parameters are determined based on the standard IAU 1976 
(precession) and IAU 1980 (nutation) theories, and are expressed as a function of TDT. The 
nutation in obliquity and longitude are corrected using series provided by the International 
Earth Rotation Service (IERS [15]), which determines these angles to high precision via 
regular Very Long Baseline Interferometry (VLBI) observations of distant quasars. The po- 
lar motion with respect to the International Terrestrial Reference Frame (ITRF), and small 
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variations in the length of day (i.e. , UTl— UTC) are also taken into account. Coordinates 
of earthbound DSN antennae, referred to the ITRF, are also known to centimeter precision 
or better based on VLBI, are taken from an existing DSN publication [18] (but see Sec. V). 


A. Equations of Motion 


The trajectory of the Pioneer 10 spacecraft must be determined by integrating the equa- 
tions of motion over the time interval of interest, given a trial set of initial conditions. The 
equations of motion I used were 

dv/dt = ajv + as + ap 

dv/dt = v (4) 


where a^ is due to Newtonian gravity, a 5 is the acceleration due to solar radiation pressure; 
and ap is an anomalous acceleration term (i.e., that which is not accounted for by known 
physics) . 

A02 considers additional terms for the acceleration which allow for alternate theories of 
gravity (their equation 3). I find that over the span of the data, these terms are always 
smaller than 3 x 10 -12 cm s -1 , and thus I neglect them for the purposes of Doppler tracking 
analysis. Other accelerations which I disregard: solar wind pressure (< 10~ 13 cm s -2 ); 
collisions with interplanetary dust (< 10 -12 cm s -2 , to heliocentric radii of 60 A.U. [19, 20]); 
and the gravitational attraction of the Kuiper belt (< 3 x 10~ 10 cm s“ 2 [2]). 

The Newtonian gravitational acceleration was computed as 


_ V" GMjijCj - r) 


a N = ^2 


(5) 


where Mj and r ; - are the mass and position of solar system body j, referred to the J2000 
coordinate frame. The bodies included in the sum were the Sun, moon, and planets, the 
positions of which were interpolated from the JPL DE405 planetary ephemeris [8, 22]. 

In a manner similar to A02, I model the acceleration due to solar radiation pressure as 
radially directed outward from the Sun with a magnitude 


as — 


cM P 


1 A.U. 

r-r© 


COsO 


( 6 ) 
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where r 0 is the barycentric position of the Sun, and the other constants in the equation are 
defined in Table I. For the Pioneer 10 geometric area I have used the area of the high gain 
antenna, which has a radius of 137 cm. The angle of the antenna to the Sun, 9 , is always 
less than 1.5° for Pioneer 10 after 1987, and here I have approximated it as 9 = 0° with a 
loss in precision in acceleration of < 4 x 10 -12 cm s -1 . 

The Pioneer 10 anomalous acceleration, a P , is modeled primarily as a constant accel- 
eration, ap = apf, where here r is a unit vector pointing from the Sun to the spacecraft. 
As noted by A02 (and below), the Doppler tracking data for Pioneer 10 do not permit one 
to distinguish between a geocentric or heliocentric acceleration, so this representation is 
also equivalent to an acceleration directed along the Earth-spacecraft line. As I am using 
the “usual” sign convention for frequencies and velocities [17], a negative value for ap will 
represent an apparent acceleration towards the Sun. 

I also test the constancy of the acceleration by adding a jerk term, 

a P (t) = (a P (0) +j P t)r (7) 

where jp is the anomalous jerk, which measures the deviation of the acceleration from a 
constant. This expression can be rewritten as 

a p(t) = ap(0)(l + t/Tj p ) (8) 

where Tj p = jp/ap{ 0) represents the timescale over which the anomalous acceleration 
changes. However, since the heliocentric spacecraft velocity is nearly constant with time 
(heliocentric radial velocity range of 13.1-12.6 km, with a mean of 12.8 km s -1 ), the jerk 
term is also equivalent to a spatial gradient of the anomalous force, and equation 8 can also 
be rewritten as 

a P (t) = a P (0)(1 + r/R jp ) (9) 

where r is the heliocentric distance and Rj p — jp/(ap(0)v r ) is the physical distance scale 
for variations in the acceleration. 

B. Spacecraft Maneuvers 

The Pioneer 10 antenna is designed to point towards the Earth. As the spacecraft moves 
outward through the solar system, regular maneuvers must be made to adjust the spacecraft 
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attitude to maintain an Earth-pointing direction. The spacecraft has two thruster assemblies 
mounted on the rim of the high gain antenna, which are aligned with the antenna and 
spin axes. During the maneuvers, the thrusters execute several small pulses, with each 
thruster assembly firing in opposite directions. The spin axis is gradually precessed until 
a spacecraft feedback loop determines that the antenna axis is again pointed towards the 
Earth. According to A02, the maneuver duration is about 15 minutes. 

In principle, the impulses from the thrusters are in opposite directions, and thus should 
impart no net change in velocity to the spacecraft. In practice, the control of the thruster 
nozzles is imperfect, and it is possible that a small velocity change will be imparted during 
the maneuver. In most cases, these velocity increments or decrements are directly visible 
in the Doppler tracking data (see Sec. V). I treated these velocity changes as adjustable 
parameters. For the j th maneuver, I modeled the velocity change as Avj = where 

A Vj is a free parameter and r ; - is a unit vector which points from the Earth to the spacecraft 
at the time of the maneuver. 

The precise epochs of the maneuvers are not easily determined from the ATDF data 
available from the NSSDC archive. In principle these data should always be available, in a 
“high Doppler rate” mode, since maneuvers can only be performed during tracking passes. 
Unfortunately very little of the high rate data is present in the archive. Rather than guess 
at the maneuver epochs, I requested and obtained from the Anderson group a file which 
contained the epochs of the maneuvers as used in the CHASMP program [21]. However, the 
velocity increments and directions were determined by my own independent analysis. 

Since the maneuvers are modeled with a single quantity, they determine the mean velocity 
shift per maneuver. Shorter time scale effects, like transient leakage from the thruster 
nozzles, will not be modeled. However, A02 found that transient effects were small, and I 
will not model them further. A02 provides a sample case of a maneuver from December 23, 
1993. This data was also present in the NSSDC archive files, and I was able to verify that 
the behavior was very close to that described by A02. 

C. Spacecraft Spin 

The downlinked tracking signal is affected by the spacecraft spin (equation 1). The 
nominal spin period is approximately 4.4 rotations per minute (rpm), however the actual 
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spin period has varied between between 4.25 and 4.55 rpm over the time span of the data 
considered in this paper. Like the maneuver data, it could be possible to determine the 
spacecraft spin from high rate Doppler data taken during precession maneuvers. Because 
this data was largely unavailable from the NSSDC archive, I also obtained a file from the 
Anderson group which contained a detailed Pioneer 10 spin period history [21]. These spin 
rate measurements come from a variety of sources [2], including the star sensor, the Imaging 
Photo Polarimeter, and the Doppler signal from precession maneuvers. During the analysis, 
I performed linear interpolation between tabulated points in the file. 

While I did not determine the spin history independently, there are several mitigating 
factors. First of all, a spin rate of 4.55 rpm corresponds to an effective frequency shift of 
~75 mHz. As I will show, a signal of this magnitude could in principle be detectable against 
the Doppler noise, but it is much smaller by a factor of ~ 40 than the signal due to the 
anomalous acceleration by 1994. Second, the more relevant quantity is how the spin rate 
change affects the Doppler frequency. For my data, the frequency shift due to the spin 
rate change is only ± 3 mHz, which is essentially undetectable. I checked these results by 
performing an analysis run where the spin rate was held fixed at its mean value, and also at 
zero, and the changes were negligible. 

D. Integration of Equations 

The equations of motion were integrated using an Adams-Bashford-Moulton predictor 
corrector algorithm, based on the DDEABM [23] routine of the SLATEC library [24] (trans- 
lated to IDL [8]). This integrator is of variable order (up to order 13) and adaptive step 
size. I adjusted the error control parameters so that frequency residuals were less than 0.1 
mHz. The initial conditions were the initial spacecraft position and velocity referred to the 
solar system barycenter. 

E. Additional Filtering 

Several additional data filtering criteria were applied, which relate to the effects of the 
Earth’s troposphere and the Sun’s corona on the Doppler signal. The Earth’s troposphere is 
known to introduce an additional signal propagation delay on the order of tens to hundreds 
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of nanoseconds. This effect is strongly dependent on the elevation angle of the spacecraft 
to the horizon. At low elevation angles, the secant effect multiplies the tropospheric delay 
by several times. In addition, there are terms in the tropospheric delay which depend both 
on the season, and atmospheric conditions at the time of the observation. These effects are 
most readily seen in ranging experiments were signal delay directly corresponds to range 
error. For Doppler tracking data, the tropospheric effect enters more subtly, as the time 
derivative of the delay, since during a single tracking pass the spacecraft’s apparent position 
will generally increase or decrease in elevation. 

Using my best-fit model (see below), I divided the residuals into intervals based on their 
elevation angle at time of reception. Figure 1 shows the root mean squared (rms) residuals 
in each interval, and demonstrates that below 15° elevation there is a strong increase. While 
A02 chose to apply an elevation-dependent weighting function which included data at low 
inclinations, but at a reduced weight, I simply excluded points for which either the received 
or the transmitted elevation angle was smaller than 15°. 

The solar corona also affects the quality of the data. During solar conjunctions, the 
trajectory of photons passes within ten solar radii of the center of the Sun. Similar to the 
troposphere, the bulk solar corona introduces a variable delay of 0-1.7 /as. The derivative 
of this variation may be imprinted on the Doppler signal. However A02 found that the net 
effects of the solar corona were small, and ultimately ignored them. I constructed a similar 
coronal model to that of A02 and also found that the the net effects of the solar corona 
on the Doppler signal were small. However the solar corona is not a uniform medium. In 
addition to the net propagation delays due to the coronal plasma, there is a general increase 
in the Doppler noise. Figure 1 also shows a plot of the rms residuals as a function of the 
photon trajectory impact parameter. Clearly the noise is enhanced for impact parameters 
within 7 x 10 12 cm (~ 0.5 A.U.), and so I also elected to exclude any trajectories which 
passed within that region. The excluded segments are shown in Figure 3, labeled as “C.” 

Finally, I found that there were several segments of data that were particularly noisy, 
and also elected to exclude those from further analysis. These segments were from 14 to 
29 September 1987, 18 January to 23 February 1992, and 13 March to 29 April 1992, and 
are shown in Figure 3, labeled as “N.” I could not find a direct explanation for why these 
particular segments were of a lower quality than the others. 

I should note that the exclusion of the segments mentioned above had a small effect on 
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the result. When, in a separate analysis, I included all of the data, the same value for the 
anomalous acceleration was reproduced to within 6%. However, because of the sensitivity 
of the least squares optimization technique to outliers, it is prudent to exclude highly noisy 
data which can significantly bias the result. 

F. Least Squares Optimization 

The Doppler data were fitted to the model iteratively using a least squares technique. 
The fitting code is based upon MINPACK-1 [25, 26], but translated to IDL [8]. The free 
parameters are: (1) the position and velocity of the spacecraft at the initial epoch; (2) an 
anomalous acceleration; (3) velocity increments A Vj due to maneuvers (a total of 18 incre- 
ments); and (4) in some cases a jerk term. The Earth station coordinates and velocities were 
also preliminarily considered to be free parameters. Upon completion of the fit, parameter 
uncertainties were estimated by adjusting the Doppler frequency uncertainties so that the 
X 2 value was equal to unity, and appropriately rescaling the parameter uncertainties derived 
from the covariance matrix of the fit. In addition to providing the parameter uncertainties, 
this method also provides an estimate of the variance of the Doppler residuals for a, given 
model. 

Since outliers can still be a problem, I gradually removed the outliers by applying a 
threshold filter. Initially the acceptance region for residuals was ±10 Hz around zero. As 
the fit steadily improved, I narrowed the acceptance region until I reached a minimum of ±60 
mHz. The distribution of residuals for the best fit is shown in Figure 2. The distribution has 
a clear sharp peak (la width of 4.2 mHz), with broad wings that extend at least to 30 mHz 
and beyond. Thus, the measured variance in the residuals will always be larger than 4.2 
mHz, and depend largely on the size of the acceptance window. I decided that a ±60 mHz 
window was a reasonable compromise between too lax and too aggressive outlier removal. 

A total of 312,116 Doppler records passed the preliminary filtering process described in 
Sec. III. I found this number of data points to be unwieldy to process simultaneously in 
core memory of a typical workstation computer, both in terms of memory consumption and 
processing time. I elected to literally decimate the data, taking only every tenth sample. 
This resulted in 31,211 raw records for the main processing runs. After application of the 
corona, tropospheric and low-noise selection criteria, a total of 23,852 or 76% of the records 
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remained. As a consistency check, I applied the same analysis to successive independent 
batches of 31,211 records drawn from the full pool of Doppler records. I found that each 
batch produced comparable results to the main batch. The distributions of parameter 
values from all ten batches were well matched by the error estimates taken from the rescaled 
covariance matrix, and therefore I have reasonable confidence that the covariance matrix 
produces appropriate statistical parameter uncertainties, even in the environment of outlier 
points. 


V. RESULTS 

In my best fit model I can confirm the signature of a constant acceleration acting on the 
Pioneer 10 spacecraft. Figure 3 shows the best fit model with and without the anomalous 
term [17]. Table II shows the best fitting anomalous acceleration value, ap, for various cases, 
as well as the inferred variance in the Doppler residuals. The best fitting value of ap for 
1987-1994 is (7.70 ±0.02) x 10 -8 cm s -2 , where the uncertainty is statistical only. It is clear 
that the Doppler residuals show an increasing trend. By the end of the data span in 1994, 
the frequency of the received Doppler signal is higher than expected by approximately +2.7 
Hz in a single round trip. 

The entry labeled “A02 Interval I” refers to the specific time intervals defined by A02. 
Interval I spans 3 January 1987 to 17 July 1990 and Interval II spans 17 July 1990 to 12 
July 1992. A02 also considers a third interval which continues up to July 1998. Because my 
data set contains a large gap from 1990.5-1991.5, and no data beyond 1994.3, I consider it 
inappropriate to quote a value of ap in Intervals II and III for direct comparison to A02. 
A02 finds anomalous accelerations in Interval I of 8.02 + 0.01 and 8.25 + 0.02 for the SIGMA 
and CHASMP techniques, respectively (xlO -8 cm s -2 ). Here I have used the weighted 
least squares values with no corona model, as these are the most comparable to my own. 
Generally, there is good agreement between the work of A02 and myself. 

The best fit case was performed without a jerk term. When a jerk is included, the fit 
improves slightly (as judged by the reduction in the rms residuals), and the anomalous 
acceleration value increases by about 5%. The small fitted jerk value demonstrates that the 
anomalous acceleration is reasonably constant over time, even when allowing the maneuver 
parameters to vary. 
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I have included two fits with simplified spacecraft spin models. The first model, “Mean 
spin,” assumes that the spacecraft spin remains constant at its mean value of 4.40 rpm. The 
second model, “No spin,” assumes that the spacecraft has no spin at all. Both cases produce 
results that are essentially indistinguishable from the best fit case, with similar values of 
a P and similar qualities of residuals. Thus, while the spin data was not independently 
determined by myself, it has little impact on the final result. The reason that the zero spin 
solution does not contain ~ 75 mHz residuals is that these residuals are essentially constant, 
and can easily be absorbed into the other free parameters, such as the initial velocity vector 
and the maneuver velocity increments. 

Taking this possibility to its logical extreme, one might surmise that the entire anomaly 
could be absorbed into the other free parameters. The next entry in Table II, “Only maneu- 
vers,” fixes ap = 0 while allowing the other parameters to vary. I found that the anomaly 
is indeed absorbed into the maneuver velocity increments, as might be expected. However 
this possibility is not likely for several reasons. First, the rms residuals are considerably 
worse. If one were to take the “best fit” case as a good fit, i.e. a reduced x 2 value of unity, 
then the “only maneuvers” case would have a reduced x 2 value of 7.2, which is very unlikely 
statistically. The residuals also show systematic trends which actually magnify the Doppler 
discontinuities across maneuver epochs (Figure 4). Also, one would have to explain how a 
set of maneuvers, whose times are irregularly spaced, could produce a steady increase in the 
velocity of the craft over 7.5 years to within one percent. Finally, the maneuver velocity 
impulses must be significantly larger in magnitude than the “best fit” case by a factor of 
~ 7.5. For these reasons I believe that the “only maneuvers” case to be extremely unlikely. 

The final entry in Table II is the case where no maneuvers are modeled, i.e. all of the 
A Vj are set to zero. Of course, Doppler discontinuities are clearly visible in this case (Figure 
5), but they are (a) small compared to the anomaly, and (b) both positive and negative 
sign, compared to the anomaly which is unidirectional. This result shows that even without 
any modeling of spacecraft maneuvers, the anomaly is significantly detected, and although 
considerable effort was put into accurate maneuver modeling, even a crude model would 
have sufficed. 

In all of the values cited in Table II, the positions and motions of the Earth stations 
were fixed to the values determined from VLBI [18]. [The station motions are due primarily 
to tectonic drift.] In another fit (not shown), I allowed the station coordinates to be free 
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parameters. I found that the fitted station coordinates converged to the quoted positions to 
within a few meters. Therefore, I left the stations fixed to their fiducial positions. 

While A02 divided their data set into three separate intervals, I do not believe this 
approach to be appropriate for the abbreviated data set that I have access to. Therefore, in 
my discussion I quote the “best fit” value, which covers the entire 1987-1994 range. 

A02 discovered annual and diurnal signatures in their residuals, which had amplitudes of 
approximately 10 mHz each. While the source was ultimately undetermined, A02 believed 
the periodic residuals to be due to previously unrealized errors in the tabulated solar system 
ephemerides, and therefore considered it to be a systematic uncertainty in the analysis. As 
can be seen in the bottom panel of Figure 3, I also detect modulations of the ^annual 
residuals at a similar amplitude. I also consider this effect to be a systematic uncertainty. 

Finally, I considered the geometric origin of the anomalous acceleration. As I have already 
mentioned, I assumed that the anomalous acceleration was directed toward the Sun. In a 
separate fit, I adjusted the equations of motion so that the acceleration was directed toward 
the instantaneous position of the Earth instead of the Sun. This change altered the Doppler 
residuals systematically by less than 0.5 mHz, and altered the best fit anomalous acceleration 
value by less than 2 %. Thus, the center of acceleration could be either the Sun or the Earth 
and still be consistent with the data. [27]. 

VI. ANOMALOUS ACCELERATION AND UNCERTAINTIES 

A02 presented a comprehensive discussion of the systematic uncertainties associated with 
the determination of the anomalous Pioneer acceleration. I do not intend to repeat such 
a discussion, but instead will summarize and adjust it. A02 divided the uncertainties into 
three main categories: those generated external to the spacecraft, those generated on board 
the spacecraft, and computational uncertainties. 

A02 estimated that the uncertainties associated with effects external to the spacecraft 
were essentially negligible, with an rms contribution of ~ 0.04 x 10 -8 cm s -2 . The largest 
estimated systematic uncertainties were associated with effects generated on board the Pi- 
oneer 10 spacecraft. A02 estimated the rms contribution of these effects to be 1.27 x 10~ 8 
cm s~ 2 , which included terms for the reflected heat from the RTGs; differential emissivity of 
the RTGs; non-isotropic radiative cooling of the spacecraft; gas leakage; and other smaller 
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effects. I adopt those values here. 

The third category, computational uncertainties, were estimated to be ~ 0.35 x 10 -8 cm 
s~ 2 , and included terms for consistency of modeling (<7 cons ist-modei) and the unmodeled annual 
and diurnal residuals. A02 was able to rely on their Interval III (July 1992-July 1998) for 
the most consistent determination of anomalous acceleration, but most of that data was 
not available to me. Thus the consistency between different models in my analysis will by 
necessity be less. For the purposes of this work, I will take cr consist . mode i to be one half of the 
range of anomalous acceleration determinations, or 

•^consist-model = 0.21 X 10 Cm S , (10) 

where, to be conservative, I have included the “extreme” cases in Table II (compare to a 
value of 0.13 x 10 -8 cm s -2 determined by A02). Thus, the total estimated computational 
uncertainty is ~ 0.38 x 10 -8 cm s -2 . The combination of the uncertainties from all three 
categories, assuming they are are uncorrelated, is op = 1.34 x 10 -8 cm s -2 . 

A02 also identified experimental “biases,” which were other effects that would tend to 
systematically increase or decrease the anomalous acceleration from its experimentally de- 
termined value. For example, they estimated that the radio transmitter exerts a radiation 
force which accelerates the spacecraft at 1.10 x 10 -8 cm s -2 , directed away from the Sun. 
This acceleration would tend to increase the anomalous acceleration. Their final bias value, 
using the sign convention of this paper [17], is bp = —0.90 x 10~ 8 cm s -2 , which I also adopt. 

Clearly, the uncertainty in the determination of the anomalous Pioneer 10 acceleration is 
systematics-dominated and not statistics-dominated. Determination of the absolute jerk is 
therefore similarly dominated by systematic uncertainties. Formally, I take the upper limit 
to the absolute jerk to be 

\j P \ <o P /T = 5.7 x 10" 17 cm s~ 3 , (11) 

where T is the data time span of 7.5 years. This upper limit is a factor of ~1.5 larger 
than the value determined in Table II. A more interesting quantity is the relative jerk (e.g., 
equation 8). The Doppler tracking data alone show a reasonably linear correlation with 
time, and hence require a small relative jerk. The effects of a jerk term would be strongest 
in the 1994-1998 time range of the A02 data set, but A02 did not see the effect. I will 
therefore still consider the jerk term shown in Table II, expressed relative to ap, to be an 
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upper limit. This leads to 

\jp/ap\ = T jp > 70 yr (12) 

and 

R jp > 170 A.U., (13) 

which implies that the anomalous acceleration, if it varies, must do so on broad spatial or 
temporal scales. 

For my determination of the anomalous acceleration I will assume that the jerk is zero, 
and hence use the “best fit” case of Table II. Following the terminology of A02, 1 label that 
experimentally determined quantity to be ap( expe r ) = —7.70 x 10~ 8 cm s -2 . After adding the 
bias value and assigning the systematic uncertainty, I arrive at 

ap = Op(exper) + bp ± Up (14) 

= (-8.60 ± 1.34) x 10~ 8 cm s" 2 . (15) 

VII. DISCUSSION 

My best value of the anomalous acceleration agrees quite closely with the value deter- 
mined by A02 (ap(A 02 ) = (-8.74 ± 1.33) x 10 -8 cm s~ 2 ), although it should be pointed out 
that I have essentially adopted their error analysis estimates directly. The rms residuals of 
all of the “non-extreme” cases in Table II are of order 8 mHz. This variance level is half of 
the standard error of 15 mHz that A02 assigned to their Doppler data processing, and thus 
compares quite favorably with their result (despite the outliers). 

The scope of this paper is to verify the Pioneer 10 anomalous acceleration by perform- 
ing an independent analysis. I will however discuss briefly some implications for alternate 
explanations of the effect. 

A02 mentions the Yukawa potential [28] as a candidate form of modified gravity, 


U(r) 


GM n+ae-’^ \ 
r \ 1 + a ) 


(16) 


where a and A are adjustable parameters. Upon computing the gradient of this potential 
and expanding in a power series of heliocentric radius r, one finds 


a GM 
1 + a 2A 2 


2 r 
1 ” 3 A 
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where the first term is the Newtonian acceleration. The second term has a clear analogy to 
equation 9, where ap(0) is the term in parentheses, and the length scale A = 2Rj p /3 > 110 
A.U. Thus, if the Yukawa acceleration — or any other modified-gravity acceleration — were 
to deviate from Newtonian plus a constant, this deviation would occur over spatial scales 
larger than the planetary solar system. 

The anomalous acceleration has been proposed to be caused by radiation from the RTGs 
or electronics in the instrument compartment [29-32]. All electric power on-board is derived 
from the RTGs, which in turn derive their power from radioactive decay of 238 Pu, with a 
half-life of r = 87.74 yr. This radioactivity also produces considerable waste heat of ap- 
proximately 2000 W. As little as 63 W of electromagnetic radiation, if radiated directionally 
away from the Sun, could explain the anomalous effect. A02 has advanced several arguments 
against these classes of explanations for the anomaly. One argument is that the anomaly is 
well enough determined over time that the radioactive decay of 238 Pu should be detectable, 
but is not. If the acceleration were related to heat dissipation, then its functional form 
would be 

a P = ap(0)2 -f / T ~ Op(0)(l - tln2/r). (18) 

This equation is again a direct analog of equation 8, however with the constraint that 
t = Tj p In 2 > 50 yr. This limit still accommodates the half-life of 238 Pu, so an explanation 
based on radiation from the RTGs cannot necessarily be excluded by my analysis of the 
1987-1994 Pioneer 10 data [33]. 

If the Doppler errors are considered to be approximately constant over time, then the error 
in the jerk should scale as T -2 , so additional data over a longer baseline could and should 
be much more constraining. A02 considered the constancy of the anomalous acceleration 
by dividing the data into three separate intervals, and attempting to analyze the intervals 
independently of one another. They found a variation of 2.0-5. 1% between their Intervals 
I and III. The “jerk” solution I present here would produce a variation in the acceleration 
between the midpoints of the two intervals of about 8%, a value which is not unreasonably 
inconsistent with the results of A02. 

In order to test the sensitivity to a jerk term, I performed a test using simulated data. I 
used the best fit trajectory with jerk to construct a synthetic Doppler series, without noise, 
over an 11.5 year baseline on a regularly sampled time grid. I then fitted that series to 
a model with no jerk, but including maneuvers. I found that a reasonably good fit could 
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be found. The rms residuals were ~ 1 mHz, which is much smaller than the typical rms 
residuals of the actual best fit models. The signature of the jerk was a small parabolic curve 
in the residuals in each segment between maneuvers. Thus, I consider it possible that a jerk 
term could be present in the residuals without being readily apparent. An analysis of the 
full Doppler data set would be desirable. 


VIII. CONCLUSION 

I have confirmed by independent analysis that the Pioneer 10 anomalous acceleration 
exists in the Doppler tracking data, and is likely not to be an artifact of the software 
processing by A02. Direct comparison to A02’s SIGMA acceleration value in their Interval 
I yields agreement at better than the 1% level. The anomaly is robust to different choices 
of spacecraft spin model, and also produces a consistent value even when all maneuvers 
are removed. This data does not constrain whether anomalous acceleration is geocentric or 
heliocentric. By including a jerk term, I have showed that the acceleration is reasonably 
constant as a function of time over a 7.5 year time baseline, but not constant enough to rule 
out thermal radiation effects due to radioactive decay of Plutonium on board the spacecraft. 
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TABLE I: Adopted Pioneer 10 and Solar Parameters 


Parameter 


Value 

Pioneer 10 Mass, Mp (g) 


2.51883 x 10 5 

Pioneer 10 Area, Ap (cm 2 ) 


5.90 x 10 4 

Solar Radiation Constant, /© (erg cm -2 s -1 ) 


1.367 x 10 6 

Reflectivity Coefficient, K, 


1.77 

TABLE II: Pioneer 10 Anomalous Acceleration (Various Procedures) 

Description 

ap 

RMS residuals 0 


10 -8 cm s' 2 

mHz 

Best fit 

-7.70 + 0.02 

7.9 

A02 Interval I 

-7.98 ± 0.02 

7.1 

With jerk ( jp = (+3.7 ± 0.2) x 10~ 17 cm s~ 3 ) 

-8.13 + 0.02 

7.8 

Mean spin (/ sp i n = 4.40 rpm) 

-7.72 + 0.02 

7.9 

No spin (/ S pi n = 0) 

-7.74 ± 0.02 

7.9 

Only maneuvers (ap = 0) 

*© 

o 

o 

o 

21.3 

No maneuvers ( { Avj } = 0) 

-8.10 + 0.01 

30.2 


a Assumes a window of ±60 mHz 
^Parameter was fixed 
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FIG. 1: Plot of RMS residuals as a function of elevation of spacecraft above Earth horizon at 
time of reception (top) and as a function of point of nearest approach to the Sun of the photon 
trajectories (bottom). Data to the right of the vertical dotted line were used in the final analysis. 
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FIG. 2: Distribution of residuals of the best fit model for only the filtered data (solid line) and 
for all of the data (dashed line). The curve is characterized by a sharp central peak, well fit by a 
Gaussian with the width shown. The distribution contains significant tails. 
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FIG. 3: Doppler residuals as a function of time of the best fit model. The top panel shows the 
residuals after setting ap — 0, and demonstrates the linear increase with time. The top panel 
shows all of the data, including segments that were filtered out because of interference due to 
the solar corona (designated by a horizontal bar with “C”) or due to general noise (designated 
“N”; see text). The bottom panel shows the filtered residuals, including the best fit value of the 
anomalous acceleration. The equivalent spacecraft velocity is also shown. Velocities and frequencies 
are plotted using the “usual” sign convention [17]. 
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FIG. 4: Doppler residuals as a function of time of the “only maneuvers” case, showing how the 
maneuver parameters can absorb some but not all of the anomaly when ap is set to zero. Note the 
change in vertical scale from the bottom panel of Figure 3. 
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